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ABSTRACT 

We present a new method to calculate fundamental Doppler measurement limits with 
a dispersed fixed-delay interferometer (DFDI) in the near infrared (NIR) wavelength 
region for searching for exoplanets around M dwarfs in the coming decade. It is based 
on calculating the Q factor, a measure of flux- normalized Doppler sensitivity in the 
fringing spectra created with DFDI. We calculate the Q factor as a function of spectral 
resolution stellar projected rotational velocity V sm.i, stellar effective temperature 
Tcfj and optical path difference (OPD) of the interferometer. We also compare the 
DFDI Q factor to that for the popular cross-dispersed echelle spectrograph method 
(the direct echelle (DE) method). We find that: 1, Qdfdi is a factor of 1.5 to 3 
higher than Qde at R ranging from 5,000 to 20,000; 2, Qdfdi and Qde converge at 
a very high R {R >100,000); 3, Qdfdi increases as R increases and Vsini deceases; 
4, for a given R, Qdfdi increases as T^s drops from 3100K to 2400K (M4V to M9V). 
We also investigate how Qdfdi is affected by OPD and find that a 5 mm deviation 
from the optimal OPD does not significantly affect Qdfdi (10% or less) for a wide 
range of R. Given the NIR Doppler measurement is likely to be detector-limited for a 
while, we introduce new merit functions, which is directly related to photon-limited RV 
uncertainty, to evaluate Doppler performance with the DFDI and DE methods. We find 
that DFDI has strength in wavelength coverage and multi-object capability over the DE 
for a limited detector resource. We simulate the performance of the InfraRed Exoplanet 
Tracker (IRET) based on the DFDI design, being considered for the next generation 
IR Doppler measurements. The predicted photon-limited RV uncertainty suggests that 
IRET is capable of detecting Earth-like exoplanets in habitable zone around nearby 
bright M dwarfs if they exist. A new method is developed to quantitatively estimate 
the influence of telluric lines on RV uncertainty. Our study shows that photon-limited 
RV uncertainty can be reached if 99% of the strength of telluric lines can be removed 
from the measured stellar spectra. At low to moderate levels of telluric line strength 
removal (50% to 90%), the optimal RV uncertainty is typically a factor of 2-3 times 
larger than photon-limited RV uncertainty. 
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Introduction 



There are over 500 detected exoplanets as of June 2011 and about 80% of them were first 
discovered by the radial velocity (RV) technique using ground-based Doppler instrument^. The 
popular Doppler instruments are based on the cross-dispersed echelle spectrogaph design, which 
we called the direct echelle (DE) method. In this method, the RV signals are extracted by di- 
rectly measuring the centroid shift of stellar absorption lines. The fundamental photon-limited 
RV un certainty usi n g the DE method ha s been studied and reported by several research groups 



(e.g., iButler et al.l (|l996l ): iBouchy et al.l (|200ll )). The best achieved RV precision with the DE 



instrument is less than 1 m • RV precision with bright and stable stars using the High Accu- 
racy Radi al Velocity Planet S earcher (HARPS) on the European Southern Observatory 3.6 meter 
telescope (|Bouchv et al.l 12003 ). While DE is the most widely adopted method in precision RV 
measurements, a totally different RV method using a disperse d fixed delay interferometer (DFDI) 
has also demon strated its capability in discovering exoplanets (|Ge et al.ll2006l : [Fleming et al.ll2010l : 
Lee et al.ll201ll ). In this method, the RV signals are derived from phase shift of the interference 



fringes created by passing stellar absorption spectra through a Michelson t ype interferorneter with 
fixed o pt ical path diffe rence (OPD) between the two interferometer arms (jErskine &: Gell2000l : iGe 



2002a 



bl: lErskindl2003l ). The stellar fringes are separated by a post-disperser, which is typically a 



medium-resolution spectrograph. Doppler sensitivity of DFDI can be optimized by carefully choos- 
ing the optical path difference of the interferometer. The D FDI meth od is promising for i ts low 
cost, compact size and potential for multi-object capability (|Gd l2002bl ) . Ivan Evken et al.l (|2010l ) 
discussed the theory and application of DFDI in detail. However, its fundamental limit for Doppler 
measurements has not been well studied before. In this paper, we introduce a method to calculate 
photon-noise limited Doppler measurement uncertainty in the near infrared (NIR) wavelength re- 
gion, where we plan to apply the DFDI method for launching a Doppler planet survey around M 
dwarfs. 

IR Doppler planet surveys are very important to address planet characteristics around low mass 
stars, especially M dwarfs. M dwarfs emit most of their photons in the NIR region. Due to lack 
of NIR Doppler techniques, only a few hu ndreds of the brightest M dwarfs hav e been searched for 



exop l anets using optica l DE instruments (jWright et al 



200S 



2004 



Endl et al. 



20061 : 



Zechmeister et al 



Blake et al.ll2010l ). To date, only about 20 exoplanets around M dwarfs have been discovered 



compared to more than 500 exoplanets discovered around solar type stars (i.e., FGK stars) despite 
of the fact that M dwarfs account for 70% stars in local universe. Nonetheless, searching for planets 
around M dwarfs is essential to answer questions such as the dependence of planetary properties on 
the spectral type of host stars. In addition, the smaller stellar mass of M dwarfs favors detection 
of rocky planets in habitable zone (HZ) using the RV technique. However, t he stellar absorption 
lines in NIR are not as sharp as those in the visible band. Recent study by iReiners et al.l (120101 ) 
shows that precision RV measurements can only reach better Doppler precision in the NIR than in 



^http://exoplanet.eu/; http:/ /exoplanets. org/ 
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visible wavelength for M dwarfs with stellar types later than M4. In this paper, we report results 
from our study on fundamental limits with the NIR Doppler technique using the DFDI method. 

This paper is organized as follows. In ^ we briefly review the principle of DFDI. In ^ 
we describe the methodology of simulations in which we compute photon-limited RV uncertainty 
for both DE and DFDI. In ^ we compare the RV uncertainties of DE and DFDI under various 
conditions. In ^ we predict the photon-limited performance of a DFDI-based Doppler instrument 
in the NIR and discuss its potential contributions to exoplanets search around M dwarfs. In ^ 
we discuss the influence of telluric lines on precision RV measurements and different methods of 
removing telluric line effects. A summary of the study is given in ^ 



Brief Review of DFDI 



The theory of DFDI has been discussed by several papers (|Gell2002bl : lErskinelboOsI : | Ivan Evken et al 



2010 ). We briefly introduce the principle of DFDI in this section, readers may refer to previous 
references for more detailed discussion. DFDI is realized by coupling a fixed delay interferometer 
with a post-disperser (Fig. [TJ. The resulting fringing spectrum is recorded on a 2-D detector. 
The formation of the final fringing spectrum is illustrated in Fig. [2l B[v,y) is a mathematical 
representation of the final image formed at the 2-D detector and it is described by the following 
equation: 



hv 



xIT{iy,y) 



LSF{v, R), 



(1) 



where 5*0(1^) is the intrinsic stellar spectrum and v is optical frequency. Sq is divided by hu to 
convert energy flux into photon flux. IT is the intensity transmission function (Equation ([2])), y 
is the coordinate along the slit direction which is transverse to dispersion direction, (8) represents 
convolution and LSF is the line spread function of the post-disperser which is a function of and 
spectral resolution R. In Equation ([2]): 7 is visibility for a given frequency channel, the ratio of half 
of the peak- valley amplitude and the DC offset, which is determined by stellar flux So{i'); c is the 
speed of light; and r is the optical path difference (OPD) of the interferometer which is designed to 
be tilted along the slit direction such that several fringes are formed along each v channel (Middle, 
Fig. [2]). We assume the LSF is a gaussian function (Equation dS])), = z//i2/2.35 because we 
assume that one resolution element is equal to the FWHM of a spectral line. 



IT = 1 + 



LSF{uo,Av) 



cos 



/ 27rz/T(z/, y) 



27rzy2 



exp 



(^ - i^af 
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(2) 



(3) 



Fig. [3 shows high-resolution (0.005 A spacing ) synthetic spectra of M dwarfs with solar 
metalhcitv (jHauschildt et al.lll999l : lAUard et allboOlh . Teg ranges from 2400K to 3100K, and log 5 



is 4.5. No rotational broadening is added in the spectrum. Most absorption lines are shallow 
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with FWHMs of several tenths of an A. Since RV information is embedded in the slope of an 
absorption line, sharp and deep lines contain more RV information than broad and shallow lines. 
Mathematically, the slope is the derivative of flux as a function of optical frequency, i.e., dSo/dv. 
The power spectrum of dSo/dv is obtained by Fourier transform. According to properties of 
Fourier transform, J^ldSQ/du] = (ip) ■ J-'[Sq], where manifests Fourier transform, i is the unit 
of imaginary number and p is the representation oi v/c in Fourier space. We plot J-{[dSQ/dv] in 
Fig. |4] as well as the spectral response function (SRF), which is J-[LSF]. SRF at i? = 5,000 
drops drastically toward high spatial frequency (high p value) such that it misses most of the RV 
information contained in stellar spectrum. As R increases, SRF gradually increases toward high 
p where the bulk of RV information is stored. A spectrograph with R of ~100,000 is capable of 
nearly completely extracting RV information. Unlike DE, DF DI can shift J-\dSQ/dv] by an amount 
determined by the OPD of the interferometer (jErskindl2003l ). For example. Fig. Ualso shows the 
power spectrum of [dSo/dz/] of a fringing spectrum obtained with a DFDI instrument with a 20 
mm optical delay, which shifts J-[dSQ/v] by 20 mm. In this case, RV information has been shifted 
from the original high spatial frequencies to low spatial frequencies which can be resolved by a 
spectrograph with a low or medium R in DFDI. 



3. Simulation Methodology 



In the DE method , an efficient way based on a spectral quality factor {Q) was introduced 
by iBouchv et al.l (j200ll ) to calculate the fundamental uncertainty in the Doppler measurements. 
The Q factor is a measure of spectral profile information within a given wavelength region considered 
for Doppler measurements. Here we develop a similar method to calculate Q values for the DFDI 
method. Instead of representing the spetral line profile information in the DE method, the Q 
factor in our DFDI method represents stellar fringe profile information. We use high resolution 



(0.005 A spacing) synthetic stellar spectra generated by PHOENIX codel Hauschildt et al 



Allard et al. 



2001 



1999 



because observed spec tra of low mass stars do not have high enough resolution 
and broad effective temperature coverage. iReiners et al.l (|2010l ) have conducted several comparisons 
between synthetic spectra generated by PHOENIX and the observed spectra. They concluded that 
the synthetic spectra are accurate enough for RV measurement uncertainty calculation. We used 
synthetic stellar spectra of solar abundance with Tefj ranging from 2400K to 3100K (corresponding 
spectral type from M9V to M4V) and a surface gravity log (7 of 4.5. The Q factor is calculated for a 
series of 10 nm spectral slices from 800 nm to 1350 nm. We artificially broaden spectra with V sin i 
from km • s~^ to 10 km ■ s~^ assuming a limb darkening index of 0.6, which is a typical value 
for an M dwarf. We convolve the rotational broadening profile with each spectral slice of 10 nm 
to obtain a rotationally-broadened spectrum. We assume a Gaussian LSF which is determined by 
spectral resolution R (Equation ([3])). After artificial rotational broadening and LSF convolution, 
we rebin each spec tral slice ac c ordin g to 4.2 pixels per resolution element (according to the optical 
design of IRET by lZhao et al.l (I2OIOI )) to generate the final 2D image on a detector based on which 
we compute the Q factor. 
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3.1. Photon-limited RV Uncertainty of DE 



Bouchv et al.l (j200ll ) described a method of calculating the Q factor for the DE method. 



We briefly introduce the method here and the reader can refer to iBouchv et alj (j200ll ) for more 
details. Let So{i') designate an intrinsic stellar spectrum. Aq, a digitalized and calibrated spectrum, 
is considered as a noise-free template spectrum for differential RV measurement, which is related 
to 5o(i^) via the following equation: 

So{u 



hv 



LSF{jy) 



(4) 



where i is pixel number and So is divided by hi' to convert energy flux into photon flux. Another 
spectrum A is taken at a different time with a tiny Doppler shift, which is small relative to the 
typical line width of an intrinsic stellar absorption. Assuming that the two spectra have the same 
continuum level, Doppler shift is given by: 

^ = ^, (5) 

c u 

where c is speed of light and u is optical fre quency. The overall RV uncertainty for the entire 



spectral range is given by (jBouchv et al 



2001 



SVr 



Q- 



-1/2 



(6) 



where Q is defined as: 



W{t) 



A. 



oW 



1/2 



(7) 



and W{i) is expressed as: 



W{i) 



(8) 



The Q factor is independent of photon flux and represents extractable Doppler information given 
an intrinsic stellar spectrum and instrument spectral resolution R. According to Equation ([6]), we 



can calculate photon-limited RV uncertainty given the Q factor and photon flux N^. 
within the spectral range. 



3.2. Photon-limited RV Uncertainty of DFDI 

A new method of calculating the Q factor for DFDI is developed and discussed here. After 
a digitalization process, a 2-D flux distribution expressed by Equation ([T|) is recorded on a 2- 
D detector in DFDI. The digitalization process involves distributing photon flux into each pixel 
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according to: 1) pixels per resolution element (RE); 2) spectral resolution; 3) number of fringes 
along slit. BQ{i,j), which is a noise-free template, is then calculated. B{i,j) is a frame taken at a 
different time with a tiny Doppler shift, i is the pixel number along the dispersion direction, and j 
is the pixel number along the slit direction. The observable intensity change at a given pixel (i, j) 
in DFDI is expressed by: 



B{i,j)-Bo{i,j) 



du{i) 



5v{i) 



dBo{i,j) 6v 

... z^U)- 

ai/(?) c 



(9) 



The Doppler shift is measured by monitoring the intensity change at a given pixel in the equation: 

6v B{i,j)-Bo{i,j) 



dBo{i,j) 
du(i) 



(10) 



Frame Bq is assumed to be a noise-free template and the noise of frame B is the quadratic sum of 
the photon noise and the detector noise ao- 



(11) 



Equation (jlip is approximated under photon-limited conditions as Brmsihj) = \/B{i,j). There- 
fore, the RV uncertainty at pixel is given by: 



(12) 



The overall RV uncertainty for the entire spectral range is given by: 



SVr 



-1/2 



-1/2 



1 



-1/2 



(13) 



where 



W{i,j) = 



B{i,j) ■ 



(14) 
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and 



nl/2 



Q 



(15) 



Equation (|15|) calculates the Q factor for the DFDI method, which is also independent of flux and 
represents the Doppler information that can be extracted with the DFDI method. According to 
Equation (jl3p . we can calculate photon- limited RV uncertainty given the Q factor and photon flux 
A^e- = / , BQ{i, j) within the spectral range. 



Optical Path Difference (OPD) of a fixed delay interferometer is a crucial parameter that affects 
the Doppler sensitivity of a DFDI instrument (Ge 2002). An optimized OPD can help increase 
the instrument Doppler sensitivity. We calculate the optimal OPD for spectra of various Tpfj and 
V sin i at different spectral resolutions (Table [T]) . We assume a wavelength range from 800 nm to 
1350 nm and an OPD range from 10mm to 41mm with a step size of 1mm in the calculation as 
described in ^3.2i Optimal OPD is selected as the one which results in the highest Q factor value. 
Increasing V sini or decreasing R naturally broadens absorption lines, decreasing the coherence 
length of each stellar absorption line (Ge 2002). Consequently, our simulations show in general 
that the optimal OPD decreases with increasing Vsini or decreasing R values (Fig. [5]). We also 
note that T^s influence on optimal OPD is not significant. 

We investigate how the Q factor is affected if OPD is deviated from the optimal value. We 
calculate the Q factor when the actual OPD is deviated from the optimal OPD by 5mm. We choose 
the lower value of the two Qs from the 5 mm deviation from the optimal delay (both positive and 
negative sides) as Qdeviated- We plot the ratio of Qoptimai and Qdeviated as a function of spectral 
resolution R in Fig. [6l We found that deviating OPD by 5mm does not result in severe degradation 
of the Q factor. The maximum degradation is 1.115 and occurs at R of 5,000 and V sin i of 5 km • s~^ 
for a star with T^s of 2800 K. The degradation can be compensated by increasing the integration 
time by 24% (1.115^) to reach the same photon-limited Doppler precision according to equation 
(13). The degradation becomes smaller as R increases. As shown in Fig. HI DFDI shifts the power 
spectrum of dSo/dv by an amount determined by the interferometer OPD so that the SRF has a 
reasonable response at a region where most of the RV information is contained. The SRF broadens 
as R increases. Therefore, it can still recover most of the RV information in a stellar spectrum 
even if OPD is deviated from the optimal value. At low and medium resolutions (5,000 to 20,000), 
Q factor degradation becomes larger as V svni increases. This is because rotational broadening 
removes the high frequency signal from a stellar spectrum, which makes the region containing most 




4. Comparison between DE and Optimized DFDI 



4.1. Optimized DFDI 
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of the RV information more sensitive to the choice of OPD as the power spectrum distribution 
becomes narrower due loss of high p component. 



4.2. Influence of Spectral Resolution R 

In theory, a spectrograph with an infinitely high resolution would be able to extract all the 
RV information contained in a stellar spectrum. However, in pratice, it is impossible to completely 
recover the RV information with a spectrograph with a finite spectral resolution whose spectral 
response function drops at the high spatial frequency end. Although the power spectrum of the 
derivative of the stellar spectrum is shifted to the low frequency region where most of the RV 
information is carried, the power spectrum is still broad in the spatial frequency (p) domain (see 
Fig. H]). Thus, high R can help to extract more RV information. In a wavelength coverage from 800 
nm to 1350 nm, we calculate Q values for stellar spectra with V sm.i of ,2 ,5 and 10 km • s~^ at 
different R (5,000 to 150,000 with a step of 5,000) in order to investigate the dependence of Q on i? 
(Fig. [7|). We find that more RV information (higher Q factor) can be extracted as R increases. Q 
factors for DFDI and DE converge at high R because the spectral response function is wide enough 
in the p domain to cover the region rich in RV information, not affected by the power spectrum 
shifting involved in DFDI. In addition, the Q factor at a given R increases as T^s drops from 3100K 
to 2400K, which is largely due to stronger molecular absorption features in the I, Y and J bands 
(see Fig. ED. 

We divide R into three regions, low resolution (5,000 to 20,000), medium resolution (20,000 
to 50,000) and high resolution (50,000 to 150,000). We use a power law to fit Q for both DFDI 
and DE as a function of R. The power indices x of three regions for T^s = 2400i^r are presented in 
Table [2j At low R region, x remains roughly a constant for km • s~^ < V smi < 5 km • s^^, but 
it drops for stars with V smi of 10 km • s~^ indicating stellar absorption lines begin to be resolved 
even at low R. At higher R regions, x decreases as V sm.i increases, a reduced value of x implies 
diminishing benefit brought by increasing R. Stellar absorption lines are broadened by stellar 
rotation, and they are resolved at a certain R beyond which increasing R does not significantly 
gain Doppler sensitivity. Overall, x for DE is larger than that of DFDI, especially for low and 
medium R. In other words, Qdfdi is less sensitive to a change of i?, and the DFDI instrument 
can extract relatively more Doppler information at low or medium spectral resolution than the DE 
method. For example, for slow rotators (ysini=2 km • s~^) at the low R region (R=5, 000-20, 000), 
Qdfdi oc R^'^^ . Doppler sensitivity 5vrm.s is inversely proportional to two factors: Q and y/N^.- 



according to Equation ([6]) and (jlSp . where A'^g- is the total photon count collected by the CCD 
detector. N^- cx {S/N)"^ • A^pixei, where S/N is the average signal to noise ratio per pixel, and A^pixci 
is total number of pixels. Note that N^.- (x R \i the wavelength coverage, S/N per pixel and the 
resolution sampling are fixed. Therefore, Svrms oc i^~0-63-o.5 _ ^-i i3 fQ,, in comparison, 

Svrms oc R~^'^'^ for DE given the same wavelengt h coverage and S /N per pixel. The power law is 



consistent with the previous theoretical work by iGd (|2002bl ) and lErskind (|2003l ). 
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We compare Q factors for both DFDI and DE at given R values, and the results are shown in 
Fig. [8j For very slow rotators (0 km • s^^ < Vsini < 2 km • s~^), the advantage of DFDI over DE 
is obvious at low and medium R (5,000 to 20,000) because the center of the power spectrum of the 
derivative of the stellar spectrum is at a high frequency domain which cannot be covered in DE 
due to the limited frequency response range of its SRF at low and medium R. The improvement of 
DFDI is ~3.1 times (i?=5,000), ~2.4 times (ii=10,000) and ~1.7 times (ii=20,000) respectively. 
In other words, optimized DFDI with R of 5,000, 10,000 and 20,000 is equivalent to DE with R 
of 16,000, 24,000 and 34,000 respectively in terms of Doppler sensitivity for the same wavelength 
coverage, S/N per pixel and spectral sampling (otherwise, see more discussions in H4:.'6\ the gain 
with the DFDI would be more significant for a fixed detector size and exposure time). Overall, DE 
with the same spectral resolution as DFDI at R=5, 000-20, 000 requires ~3-9 times longer exposure 
time to reach the same Doppler sensitivity as DFDI if both instruments have the same wavelength 
coverage and same detection efficiency (i.e., N^- is the same). The improvement of DFDI at 
i?=20, 000-50, 000 is not as noticeable as at the low R range. The difference between DFDI and 
DE becomes negligible when R is over 100,000. In other words, the advantage of the DFDI over 
DE gradually disappears as R reaches high resolution domain (R > 50,000). In addition, the 
improvement for relatively faster rotators (5 km • < V sin i < 10 km • s~^) with DFDI is less 
significant than it is for slow rotators. 



4.3. Influence of Detector Pixel Numbers 

In the NIR, detector pixel number is typically smaller than the optical detector. Furthermore, 
the total cost for an NIR array is much higher than an optical detector with the same pixel number. 
In the foreseeable future, detector size may be one of the major limitations for Doppler sensitivity 
improvement. We study the impact of the limited detector resource on the Doppler measurement 
sensitivity. Using the same detector resource, we find that it is fair to compare their Doppler 
performance for the same target with the same exposure to understand strength and weakness 
for each method although DFDI and DE are totally different Doppler techniques. According to 
Equation ([6]) and ()13p . we define a new merit function, 



Q' = Q- VK^, (16) 

to study photon-limited Doppler performance for both methods with the same detector size. Note 
that the newly defined merit function is directly related to photon-limited RV uncertainty, i.e., 
inverse proportionality. N^- is calculated by Equation (fT7|l : 



2.512™-^ ' 

in which F^, is the photon flux in the wavelength coverage region AA of an rrij = star with the 
unit of photons • • cm~^; rj is instrument total throughput; Stei is the effective surface area of 
the telescope; ioxp is the time of exposure; and nij is the J band magnitude. 
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Here we use IRET as an example to illustrate strengths of the DFDI method for Doppler 
measurements. IRET adopts the DFDI method and has a wavelength coverage from 800 nm to 
1350 nm and a spectral resolution of 22,000. For a fixed detector size (i.e., total number of detector 
pixels) and fixed number of pixels to sample each resolution element, the total wavelength coverage 
of a Doppler instrument, A A, is 

AA = ^^-^^, (18) 

border R-Ns' ^ ^ 

where A'pix is total number of pixels available on a CCD detector and Ns is the number of pixels 
per resolution element, Ac is the central wavelength and -Porder is the number of pixels sampling 
each pixel channel between spectral orderj^. Equation (jl8p shows that AA is inversely proportional 
to R. Table [3] gives the relation of R and A A assuming A^'pixci, -Pordcr said Ns as constants. Ac is 
set to be 1000 nm because it is approximately the center of the Y band. We calculated the ratio of 
Qdfdi ^^'^ Q'de which we use the photon flux of a star with a T^s of 2400K (Fig. [9|). Qiret 
consistently higher than Qde regardless of R of the DE instrument. In other words, IRET is able to 
achieve lower photon-limited RV uncertainty compared to a DE instrument with the same detector. 
The result seems to be different from the conclusion we drew in ^4.21 in which we compare Q factors 
of the same R and AA and reached a conclusion that DFDI with an R of 22,000 is equivalent to DE 
with an R of 35,000 (a factor of 1.6 gain) for the same wavelength coverage and the same detection 
efficiency (Fig. [7]). The key difference between this case and the earlier case is the fixed detector 
resource instead of fixed total collected photon numbers. Since lower spectral resolution allows to 
cover more wavelengths, more photons will be collected for the same instrument detection efficiency 
for both DFDI and DEM. Note that Q' consists of two components, Q and N^,- . For a given number 
of pixels on the detector, A'^g-^DFOi is higher than A^'c-^de due to the larger wavelength coverage. In 
addition, Qdfdi(AAdfdi) is more than (5de(AAde)- Consequently, we see in Fig. [9] that Qipiet i^ 
higher than Q'de all i? of a DE instrument. Fig. [9] also shows that the minimum of Qdfdi/^de 
is dependent of Vs'mi. The ratio of Q' reaches a minimum (Q^g reaches a maximum) around an 
R of 50,000 for slow rotators {Vsini < 5 km-s~^). It increases at the low R end because the 
spectrograph has not yet resolved stellar absorption lines. On the other hand, the ratio increases 
at the high R end because of fewer photons (see Table [3]). For fast rotators (F sini=10 km • s~^), 
the ratio reaches a minimum around R of 30,000. 

For a slow rotator (ysini=2 km-s^^) at low R region (R=5,000-20,000), Qdfdi « R^-^^. 
Since Doppler sensitivity Svrms is inversely proportional to two factors: Q and \/N^- according to 
Equation ([6|) and (jl3p . the Doppler sensitivity becomes nearly independent of spectral resolution for 
the DFDI method (oc R'^-^'"^) if the detection size (or total number of pixels) is fixed. This indicates 
that we can use quite moderate resolution spectrograph to disperse the stellar fringes produced by 



^In principle, each frequency channel for both DFDI and DE instruments can be designed identical. In practice, 
the DFDI instrument tends to use ~ 20 p ixels to sample f r inges in the slit direction, which is only for measurement 



convenience, not a requirement. In fact, iMuirhead et al.l (j2011h has demonstrated a phase-stepping method which 



does not require sample fringes in the slit direction.^ 
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the interferometer in a DFDI instrument while maintaining high Doppler sensitivity. This opens 
a major door for multi-object Doppler measurements using the DFDI method as proposed by Ge 
(2002). In comparison, the Doppler sensitivity for the DE method still strongly depends on spectral 
resolution for a fixed number of detector pixels (oc R~^'^'^), indicating that higher spectral resolution 
will offer better Doppler sensivity. 

4.4. Influence of Mult i- Object Observations 

As discussed in ^4.21 and 14. 3t the DFDI instrument can be designed to have a moderate res- 
olution spectrograph coupled with a Michelson type interferometer. Moderate spectral resolution 
allows a single order spectrum or a few order spectra to cover a broad wavelength region in the NIR 
region while keeping the Doppler sensivitiy similar to a high resolution DE design which requires 
a large detector array to cover spectra from a single target. This indicates that the DFDI method 
has much greater potential for accommodating multiple targets on the same detector as proposed 
by Ge (2002) than the DE instrument. In order to evaluate the potential impact of multi-object 
DFDI instruments, we redefine the merit function Q" as: 

Q" = Q-VK^-Kbj, (19) 

where Noj,j is the number of objects that can be monitored simultaneously, and a is the index of 
importance for multi-object observations. From the perspective of photon count and S/N, multi- 
object observations are equivalent to an increase of N^-, and thus a is 0.5. However, from an 
observational efficiency point of view, Q" should be proportional to N^i^j because the more objects 
are observed simultaneously, the quicker the survey is accomplished, and a is therefore equal to 1. 

We assume a detector that covers from 800 nm to 1350 nm at /2=100,000 so that we can use 
the Q factors obtained in §4.2[ N^- is a constant since we assume identical AA. Nobj is inversely 
proportional to the number of pixels per object which is proportional to spectral resolution R 
(Equation pB|) ). Note that we do not require Nobj to be an integer because we can, in principle, fit 
a fraction of spectrum on a detector to make full use of the detector. Q"s for both DFDI and DE are 
calculated. Fig. [TOl shows the ratio of Q" and Q'^^^oo ooo DFDI under two different assumptions 
of a. For a=0.5, i.e., increase of N^bj is equivalent to photon gain, only a slight improvement is 
achieved if the detector is used for multi-object observations at lower resolution than 100,000. In 
comparison, from a survey efficiency point of view (i.e., a=l), we see a factor of ~4-6 times boost 
of Q" in multi-object observations. The truncation at -R=5,000 is due to a practical reason that 
a lower resolution than 5,000 is rarely used in planet survey using RV techniques. On the other 
hand, similar calculation is also conducted for the DE method (Fig. [TT]l . in which we find that 
high resolution single object observation is an optimal operation mode for DE from a perspective 
of photon gain (a=0.5). At a=l, the increase of Q" is a factor of ~3 at the most. 

We compare the maximum of Q" for both DFDI and DE at different V sin i in Table HI At 
a=0.5, the advantage of DFDI over DE is ~1.1 for a wide range of V sin i. In other words, from the 
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photon gain point of view, there is no significant difference between DFDI and DE in multi-object 
RV instruments. However, from the survey efficiency point of view (a=l), we see a factor of 3 boost 
of Q" in DFDI for slow rotators {V sini <2km • s~^), suggesting 9 times faster in terms of survey 
speed. For fast rotators (i.e., ysini=10km • s~^), the boost drops to 1.78. Our study confirms 
that the DFDI r nethod has an advantage for multi-object RV measurements over the DE method 
as suggested by [Gg (j2002bl ). 



4.5. Influence of Projected Rotational Velocity V sm.i 

Projected rotational velocity V svni broadens stellar absorption lines and thus reduces the Q 
factor. We carry out simulations calculating Q factors of different V sini (0 km • s^^ < Vsini <10 
km • s~^) at i?=150,000 and various Teg. We assume a wavelength range from 800 nm to 1350 nm. 
The results are shown in Fig. [T2j The Q factor decreases as V sin i increases. It is clear that slow 
rotators would be better targets to reach higher photon-limited RV precision because the spectrum 
of a slow rotator contains more Doppler information. 



5. IRET Performance and its Survey Capability 



Here we use the IRET as an example to demonstrate the Doppler sensitivity with the DFDI 
method and its capability for an NIR survey for exoplanets around M dwarfs. IRET has a spectral 
resolution R of 22,000. It works mainly in the I, Y and J bands from 800 to 1350 nm . 4.2 pixels 



samp le one resolution element and 25 pixels sample the fringes in the slit (y) direction (jZhao et al 
20ld ^. 



We conduct a series of simulations as described in ij3.2| using synthetic spectra of M dwarfs 
and the instrument specifications of IRET ( Zhao et al.ll20ld ) to seek for an optimized OPD that 
maximizes the Q factor. Fig. [13] shows results from simulations of M dwarfs at different T^s- 
It is shown that for T^ff ranging from 2400K to 3100K, the Q factor is maximized at an OPD 
of around 18.0 mm for V svni between 3 km • s~^ and 5 km-s~^. In addition to the Q factor, 
we also need photon flux information for a particular star to calculate the photon-limited RV 
uncertainty (Equation (jl3p ). which is given by Equation (jl7p . We calculate using a synthetic 
stellar spectrum and then normalize it such that the flux at 1235 nm is 195.2 photons-s~^-A~^-cm~^. 
Instru ment throughput is as shown in Fig. [141 which is estimated based on our current optical 
design (jZhao et al.ll2010l ). Note that the estimation does not consider fiber coupling efficiency, fiber 
efficiency, seeing coupling loss and atmospheric transmission. We choose the telescope aperture size 
of 3500 mm as the design reference to calculate 5tei. The telescope obscuring structure is accounted 
for in the total detection efficiency, r]. We assume a tcxp of 30 min on an mj = 9 star in the RV 
uncertainty contours. Fig. [15] shows the fundamental photon-limited RV uncertainty for a set of 
grids of Tgff and V s\ni. The photon noise limited Doppler precision can reach ^-^^1-2 m • s~^ for 
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slow rotating cool stellar objects. There is a general trend of RV uncertainty increasing from the 
low Teff and low Vsini to the high T^s and high Vs'mi. Average S/N per pixel as a function of 
wavelength is plotted in Fig. [T6j 

As of June 2011, there were only six M dwarf exoplanets discovered in the northern hemi- 



sphere (IButler et al.l 2004: ICharbonneau et al.ll2009l : iForveille et al.ll2009l : iHoward et al.ll2010l : iHaghighipour et al 



2010|; Ijohnson et al.l |2010| ) . We compare the velocity semi-amplitude K of these exoplanets and 
RV uncertainty predicted for IRET (Table [5]). All of them would be detectable by IRET un- 
der photon-limited conditions. Therefore, IRET is a suitable instrument conducting follow-up RV 
measurements of these targets to confirm the detection and detect additional planet companions 
in these systems. We also compare the velocity semi-amplitude, Kuz, of an Earth-mass planet 
located within the habitable zone (0.05 AU away from the host star) to the IRET photon-noise 
detection limit and find that the RV uncertainty of IRET is slightly larger than Kuz- It may be 
challenging to detect a possible Earth-mass planet in the habitable zone of these planet systems. 
However, IRET is capable of detecting super-Earth exoplanets in the HZ or Earth-like planets in 
the HZ around nearby brighter M dwarfs under photon-limited condition. 

Earlier study showed that NIR RV measurements gain precision for objects late r than M4 be- 



cause of relatively higher flux and increasing stellar absorption features in the NIR (iReiners et al 



2010l ). In addition, our study in earlier sections showed that stellar rotation intrinsically broadens 
stellar absorption features and results in less Doppler information, leading to larger RV measure- 
ment uncertainty with IRET. Therefore, slow rotating mid-late type M dwarfs would be major 
science targets for a planet survey with IRET. 

We conducted a literature survey to estimate how many potential targets are suitable for the 
M dwarf survey in the northern sky using IRET. We select our potential survey targets from the 
LSPM catalog, which was compiled by Lepine (2005) through a photometric survey of the entire 
northern sky. This catalog is nearly complete for stars with proper motion higher than 0.15"yr~^. 
We apply a color cut, V-J>3, to select M type stars from the LSPM catalog. Table [6] shows the 
number of M dwarf targets available brighter than certain J band magnitude. However, not all of 
the bright M dwarfs are suitable targets. In order to reach high precision (1-3 m • s~^), we only 
choose M dwarfs with V sini less than 10 km • s~^ as major targets. However, previo us study shows 



that there appears t o be a trend that the median of V sin i increases as T^ff decreases (jjenkins et al 



20091 ). According to iJenkins et al.l (j2009l ). mid-late type M dwarf V sini distribution can be fitted 
by a power law with an index of -1.12, which indicates that 41% of mid-late M dwarfs have a 
V sini less than 2 km • s~^ and 74% of them have a V sini less than 10 km • s~^. They also pointed 
out that 43% of M dwarfs in their sample exhibit an Ha emission that indicates young and active 
objects not suitable for RV survey. Statistically, there are about 18% of nearby M dwarf targets 
that are suitable for a precision RV survey (inactive and with y sini < 2 km • s~^). Based on IRET 
photon-limited performance, it is likely that IRET can reach less than 3 m • s~^ photon-limited 
Doppler precision for over 200 nearby late M dwarfs to detect and characterize Earth-mass and 
Super-Earth mass planets. 
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6. Telluric Line Effect 



6.1. Quantitative Telluric Line Effect 



Ground-based NIR observation is prone to contamination of telluric lines. Precision RV mea- 
surements in the NIR is an extreme case that requires significant disentanglement of stellar absorp- 
tion lines and telluric lines in order to reach photon-limited precision. Even in the Y band (960 
to 1080 nm) and J band (1120 to 1320 nm), where tellur ic line contamination is less severe, there 
are still many shallow telluric lines (Fig. \T7 \ iLordl ( 19921 )) ^. Equation ([1]) should be rewritten as 
follows if telluric absorption lines are considered: 
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(20) 



where AT is the atmospheric transmission function and AA is the atmospheric absorption function. 
In Equation (|20|) . photon flux distribution on the detector, B, is comprised of a signal component 
Bs and a noise component Bj^j. Ideally, we require that the detector flux change 5B is entirely 
due to the stellar RV change 6vs- However, 6B is also partly induced by telluric line shift 6v!\f 
resulting from atmospheric behaviors. Therefore, both 6vs and 6vj\f contribute to 6B. We have two 
sets of RV measurements, 6vs + cr(0, Svrms,s) for RV of Bs and 6vn + a{0, 5vrms,N) for RV of -Bat, 
where (t(0, 5) represents random numbers following a gaussian distribution with a mean of and 
a standard deviation of 5. We weight RV measurements with the inverse square of photon-limited 
RV uncertainties. The flnal RV measurement is expressed by the following equation: 



5v 



{6vs + CJ(0, SVrms,s)) ' ^V^^^ g + {6vn + Cr{0, 6Vrms,N)) ' ^f^^s N 



5V ^ C + 5f ^ AT 
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(21) 



In practical Doppler measurements, 5vs consists two components, stellar RV and Earth's barycen- 
tric RV. We assume a constant stellar RV (i.e., no differential stellar RV). The Earth's barycentric 
motion may be predetermined and the effect is removed from the observed 2-D spectrum by shifting 
Bg according to the barycentric velocity. Therefore 6vs=^ for the Bg component. We further as- 
sume that Bn has an RV fluctuation of (t(0, 5vn,atm) because of the Earth's turbulent atmosphere. 
The measured RV uncertainty 6v is equal to: 
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(22) 



^http: / /atran. sofia.usra.edu/cgi-bin/atran/atran.cgi 



- 15 - 



Three examples are considered here to represent different level of telluric contamination to the 
precision RV measurements: 1, if 6vrms,s > Svrms^N, then 5vrms = ('^'^Ar,ATM + '^^rms.Tv)"^^^' 
and the RV measurement is dominated by atmospheric behaviors, this approximation applies 
in a wavelength region with dense telluric line distribution; 2, if 5vrms,S ^ ^Vrins,N (i-e., in 
very transparent atmosphere windows where few telluric lines exists), then 6vrms = ^Vrms,s^ 
and the RV uncertainty is limited by stellar photon noise; 3, for an intermediate situation, if 
Svrms,s c'-iid Svrms,N are identical, then we apply the same weight on both RV measurements. 

In the practical NIR spectroscopic observations, telluric contamination must be minimized in 
order to reach high Doppler precision. We investigate the effectiveness of two major ways to remove 
telluric contamination from the stellar spectra with IRET: telluric line masking and removing, and 
telluric line modeling and removing. The study results are summarized below. 



6.2. Telluric Line Masking 

Telluric absorption lines are not homogeneously distributed in the NIR spectra, instead, a large 
number of telluric lines are from relatively concentrated spectral regions. The simplest way is to 
mask and remove those severely contaminated regions while leaving the less contaminated stellar 
lines for RV measurements. This would reduce RV measurement uncertainty caused by telluric 
line contamination, leading to improved RV precision. However, if too much of the wavelength 
coverage region is masked and removed, then RV uncertainty due to photon noise would increase. 
Therefore, a balance between the uncertainty brought by telluric lines contamination and by photon 
noise must be found and an optimal RV performance can be achieved. 

We use an M9 dwarf (Tefj=2400K, ysinz=5km- s~^ and mj=9) as an example to illustrate 
how this masking technique affects the RV measurement precision with IRET. We assume a 30 
min exposure time, a wavelength coverage from 800 to 1350 nm and an instrument throughput as 
shown in Fig. [T31 We calculate photon-limited RV uncertainty dvs^rms and Sv^^rms according to 
Equation (fT3]l . We calculate Q factors, Qs and Qn, based on the two components, Bs and B]\f, in 
Equation (120p . The photon flux of Bs and B^ are calculated from the inputs of the stellar type, 
magnitude, exposure time, instrument specifications and telluric absorption properties. 

In practice, RV uncertainty of B^ is not dominated by photon-noise, instead, i t is dominated 



by atm ospheric behaviors such as wind, molecular column density change, etc. iFigueira et al 



(j2010l ) used HARPS archive data and found that O2 lines are stable to a 10 m • s~^ level over 6 
years. However, the stability of telluric lines becomes worse if we take into consideration other gas 
molecules such as H2O and CO2. We assume different levels of RV fluctuation due to atmospheric 
behaviors dvN,ATAi (5 m • s~^, 10 m • s~^ and 20 m • s~^). The uncertainty induced by atmospheric 
telluric lines is transferred to Svrms via Equation (I22p without the effort of disentangling stellar and 
telluric lines as if the uncertainty is due to photon- noise of Bs- Telluric lines with the largest slope 
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are masked at the highest priority because they cause the most RV uncertainty in the measurement. 

Fig. [18] shows the results of teUuric hne masking for an M9 dwarf under different levels of 
telluric line masking. After masking a certain portion of the most severely contaminated stellar 
spectrum, optimal RV uncertainty Vrms is 3.3 m • s^^, 3.6 m • s~^ and 3.9 m • s~^ respectively for 
S'Vn,atm of 5 m • s~^, 10 m • and 20 m • s~^. They are 2.6, 2.8 and 3.0 times worse compared to 
the fundamental photon-limited RV uncertainty of the stellar spectrum 5vrms,s^ i-e-, 1-3 m • s~^. If 
we do not apply telluric line masking at all, RV uncertainty is dominated by atmospheric behaviors, 
Svn,atm- It is expected that Svrms increases as we mask more of the stellar spectrum because less 
photon flux is available to calculate RV. Depending on values of Svn,atm, the optimal spectrum 
fraction used in RV measurement lies in between 30% and 50%. Higher value of 5vn,atm requires 
masking and removing larger fraction of stellar spectrum. It is noted that we would have expected 
Svrms to be only ~1.4 times worse than 5vrms,s after ~50% of the spectrum is masked for pure 
concern of photon flux. However, Svrms is ~3 times worse than 6vrms,s because not only N^- but 
also Q is reduced after telluric line masking. 



6.3. Telluric Lines Modeling 

In principle, the telluric lines can be modeled and removed using a telluric line spectrum tem- 
plate. This kind of template can be obtained by observing fast rotating early type stars since the 
spectrum of an early type star (early A type in particular) is nearly featureless and a g ood approxi 



matio n of a blackbody spectrum in NIR except for broad hydrogen absorption features (jVacca et al 



20031 ). However, the number of available telluric standard stars limits the application of this method. 
In addition, the spectrum of a telluric standard star is not entirely featureless, which results in some 
residuals that undermines precision RV measurements. 

As a first order estimation, we consider the simplest case in which we assume that the spectrum 
of a telluric standard star is a perfect blackbody spectrum with no spectral features which may 
otherwise be perfectly modeled and removed. Since the spectrum of a science target and the 
spectrum of a telluric line standard star are not taken simultaneously, the telluric lines' centroids 
and depths may change over time. We assume that the centroids of telluric lines remain unshifted. 
As a result, the only uncertainty is the depth of a telluric line. Equation (113p and (|20p offer insight 
in understanding the process of telluric line modeling. In the telluric line modeling, we model the 
telluric lines by observing a standard star or by forward modeling using a list of telluric lines. We 
remove telluric lines from the observed spectrum of a science target. The depths of the telluric 
lines are consequently reduced by a certain fraction depending upon the level of removal, which 
effectively reduces AA (Atmosphere Absorption). A reduced AA means less photon contribution 
in Bn and thus less uncertainty induced by atmospheric behaviors. 

Fig. [19] shows how RV uncertainty Svrms is correlated with different levels of telluric line 
modeling and removing. We assume 6vn,atm to be 10 m-s~^. Fig. [20] shows the optimal RV 
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uncertainty at different removal levels. If no removal is involved, the optimal 5vrms is 3 times worse 
than dvrms.s- At low levels of telluric line removal, i.e., a small fraction of the strength of telluric 
lines is removed, it requires to mask out about 30% to 50% of the stellar spectrum to reach an 
optimal RV uncertainty. As the removal becomes more effective, the influence caused by telluric 
lines decreases, and only a small portion of the spectrum requires masking. Therefore, the optimal 
RV uncertainty, 6vrms, reaches a smaller value. At a very high level of removal, in which 99% of the 
strength of telluric lines can be modeled and removed, i.e., the residual of telluric line modeling is 
only 1%, telluric line masking is not necessary. The entire wavelength coverage region can be used 
in the RV calculation. 6vrms approaches the fundamental photon- limited RV uncertainty, Svrms,s- 
In previous study, the residual of telluric lines has reached an rms of 0.7% using a 3 -component 



(abso rption gas, telluric lines and star) spectrum to model the observed stellar spectrum (jBean et al 



20ld ). Therefore, it is quite possible to reach photon noise limited RV precision with IRET using 



this telluric line modeling and removing method. 



7. Summary and Discussion 

7.1. Q Factors for DFDI, DE and FTS 

We develop a new method of calculating photon-limited Doppler sensitivity of an instrument 
adopting the DFDI method. We conduct a s eries of simulations based on high resoluti on synthetic 
stellar spectra generated by PHOENIX code dHauschildt et al.l[l999l : kllard et ahlboOlh . In simula- 



tions, we investigate the correlations of Q and other parameters such as OPD of the interferometer, 
spectral resolution R and stellar projected rotational velocity Vsini. We find that optimal OPD 
increases with increasing R and decreasing Vs'mi. Empirically, the optimal OPD is chosen such 
that the density of the interference combs matches with the line density of the stellar spectrum. 
Based on the simulation results, the optimal OPD is determined as the one that maximizes the 
Q factor. In fact, optimal OPDs found from empirical way and from numerical simulation are 
consistent with each other. For example, for Vsmi=0 km • s~^ and i?=50,000, simulation gives 
an optimal OPD of 30 mm. The interference comb density of an interferometer with OPD of 30 
mm is ~0.3 A at 1000 nm, which indeed matches the width of a typical absorption line after spec- 
tral blurring with R of 50,000. An independent method to calculate photon-noise limited Doppler 
measurement ui icertainty in the o ptical is being developed, and the results will be reported in a 



separate paper (| Jiang et al.ll201ll ). We have compared results from both methods and confirmed 
that both independent methods produce essentially the same results for both optical and NIR 
Doppler measurements. 

We investigate how the Q factor is affected if OPD is deviated from the optimal value and 
find that a deviated OPD (5mm) does not result in a significant Q factor degradation, which is 
mitigated as R increases. We find that the Q factor increases with R for both DFDI and DE, and 
eventually converge at very high R {R > 100,000). The convergence of DFDI and DE methods is a 
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natural consequence because the measurement method does not make a difference after the spectral 
resolution becomes extremely high. In addition, Q factors at a given R increase as T^s drops from 
3100K to 2400K, which is due to stronger molecular absorption features in NIR (see Fig. [3]). The Q 
factor decreases as V sm.i increases because stellar rotation broadens the absorption lines, leading 
to less sensitive measurement. 

We compare Q factors for both DFDI and DE at a given R. For slow rotators (0 km • s~^ < 
V sini < 2 km- s~^), DFDI is more advantageous over DE at low and medium R (5,000 to 20,000) for 
the same wavelength coverage AA. The improvement of DFDI compared to DE is ~3.1 (i?=5,000), 
~2.4 (/2=10,000) and ~1.7 (i?=20,000), respectively. In other words, optimized DFDI with R 
of 5,000, 10,000 and 20,000 are equivalent in Doppler sensitivity to DE with R of 16,000, 24,000 
and 34,000, respectively. The improvement of DFDI at R 20,000 to 50,000 is not as noticeable 
as at low R range. The difference between DFDI and DE becomes negligible when R is over 
100,000. For relatively faster rotators (5 km • s~^ < Vsini < 10 km • s~^), the improvement with 
DFDI is less obvious than it is for very slow rotators. DFDI has strength when the spectral lines 
in a stellar spectrum are not resolved by a spectrograph, which is the case for low and medium 
resolution spectrograph. Under such conditions, the fixed delay interferometer provides additional 
resolving powers for the system. After the lines are fully resolved by the spectrograph itself, the 
interferometer in the system becomes dispensable, which is the reason why we see the convergence 
of DFDI and DE at very high spectral resolution. 

Fundamental performance of a Fourier-t ransform spectr ometer (FTS) in the application of 



Doppler measurements has been discussed by iMaillardI (jl996l ). There are similarities between the 



FTS and the DFDI method, for example: 1, both methods use the interferometer as a fine spectral 
resolving element; 2, RV is measured by monitoring the temporal phase change at a fixed OPD of 
the interferometer. In DFDI method, OPD is scanned in each frequency channel because of two 
relatively tilted mirrors, and the resolution of the post-disperser in DFDI is chosen to ensure a 
reasonable fringe visibility. Therefore, the DFDI method is a extended version of the FTS method 
with a low-medium resolution post-disperser. However, one major difference between these two 
methods is that the interferometer itself is used as a spectrometer by OPD scanning in the FTS 
method while an additional spectrograph is employed in the DFDI method. The advantage of 
introducing an additional spectrograph into the system is that the visibility (or fringe contrast) is 
no longer limited by the bandpass as in the FTS case , which is the reason that the DFDI method 



can be applied in broad-band Doppler measurements. iMosser et al.l (|2003l ) discussed the possibility 
of an FTS working in broad band by introducing a low resolution post-disperser and concluded 
that the FTS method is inferior (by a factor between 1 and 2) to DE method even after employing 
an post-disperser. This conclusion should be accepted with cautions because they compared an 
FTS with a post-disperser (i?=1200) with a DE instrument with a much higher spectral resolution 
(i?=84,000), which is not necessarily a fair comparison. 

We define new merit functions (Equations (I16p and (I19p ) to objectively evaluate Doppler per- 
formance for both DFDI and DE methods. For Q', the merit function for single object observation. 
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we find that Q^pDi consistently higher than Qj^g regardless of the R of the DE instrument 
under the constraint of total number of pixels, i.e., both the DFDI and DE instrument adopt the 
same NIR detector. The DE instrument requires using a larger detector in order to reach the 
same wavelength coverage as the DFDI instrument. Note that the above conclusion is based on 
the assumption that the number of pixels per spectral order are the same for DFDI and DE. In 
practice, a DFDI instrument uses ~20 pixels to sample spatial direction, i.e., the direction trans- 
verse to dispersion direction, while ~5 pixels are usually used to sample spatial direction in a DE 
instrument. However, the ~20 pixels sampling is not a requirement for DFDI but rather for the 
convenience of data reduction. Normally, ~7 pixels sample one spatial period of a stellar fringe , 



which in principle are adequate based on a phase-stepping algorithm provided by lErskind (|2003l ). 



If the same detector is used, the spare part of the detector in DFDI can be used for multi- 
object observations. Consequently, in addition to single-object instrument, we also investigate Q", 
a merit function for multi-object RV measurement for both DFDI and DE. Different conclusions are 
reached depending on different value of a, an index of the importance of multi-object observation. 
From a pure photon gain point of view, DFDI and DE instruments have similar Q" values with 
^DFDi slightly better than QJ^^ (a factor of ~1.1). From a survey efficiency point of view, a DFDI 
multi-object instrument is 9 times faster than its counterpart using DE for slow rotating stars 
(0 km ■ s~^ < Vsini < 2 km • s~^) and ~4 times faster for fast rotators {V sini > 10 km • s~^). 



7.2. Application of DFDI 

InfraRed Exoplanet Tracker (IRET) is used as an example to demonstrate RV performance 
with the DFDI method and illustrate how the RV performance is related to other parameters. It 
is shown that for slow rotators with T^s ranging from 2400K to 3100K, the Q factor is maximized 
at an OPD of around 18.0 mm. According to predicted photon-limited RV precision, IRET is 
capable of detecting Earth-like planets in habitable zone around bright M dwarfs if they exist. 
In our simulations, we consider the photon loss due to optics, CCD detector and the telescope 
transmission loss of 20%. For the real observation, the total detection efficiency needs to include 
the atmospheric transmission loss, fiber coupling and seeing losses. It is likely that the overall 
detection efficiency is a factor of two lower, which requires 2 times longer exposure to reach the 
same Doppler precision. In addition, the real instrument may have short and long term systematics 
which affect RV measurement precision. Therefore, the performance of IRET predicted in the paper 
should be considered as an estimation in an optimistic case. 

There may be other practical concerns about the instrument using the DFDI method, most of 
them are due to the relative low spectral resolution compared to current DE instruments. First of 
all, an absolute wavelength calibration for a DFDI instrument is not as precise as a DE instrument 
with a higher spectral resolution. For example, at a spectral resolution of 22,000 for IRET, a 
li ne profile wit h a FW HM of ~0.45 A in Y band is expected. Following the method described 



m 



Butler et al.l ( 1996 ). it corresponds to 136.4 m • s ^ RV uncertainty at a S/N of 100 if only one 
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spectral line is used. iRamsev et al.l ([20101) proposed to use a U-Ne emission lamp as a wavelength 
calibration source and it has approximately ~500 lines in Y band according to their measurement. 
Therefore, after all the lines in Y band are considered, ~6 m • RV uncertainty is introduced in 
the process of absolute wavelength calibration. In comparison, a DE instrument at R of 110,000 
causes ~1.2 m- s~^ RV uncertainty in an absolute wavelength calibration. However, an absolute 
wavelength solution is only required for the DE method in order to measure RV drift due to 
instrument instability, which is measured in a different method in a DFDI instrument. It is similar to 
a stellar RV measurement, the difference is that the object is switched from a star to an wavelength 
calibration source. Vertical fringe movement of absorption or emission lines of an RV calibration 
source is measured instead of centroid movement measurement in a DE instrument. In this case, 
a DFDI instrument (e.g., IRET, ii=22,000) is equivalent to a DE instrument with R of 37,000 
in terms of Doppler measurement precision (see ^4.2p . Therefore, instrument RV drift calibration 
process introduces an RV uncertainty of ~3.5 m • s~^ for IRET in the example of a U-Ne lamp 
calibration source. In addition, the RV uncertainty can be further reduced by increasing S/N 
and number of measurement. Secondly, at a low spectral resolution, it is challenging to perform 
spectral line profile analysis and thus it requires high-resolution follow-up in order to confirm or 
exclude a possible detection. Last but not least, in a binary case in which the observed spectrum 
is blended, two approaches can be used for identification: 1, from measured RV, if the flux ratio is 
small, similar to the planet companion case, then the lower mass companion can be identified in 
the measured RV curve even though small flux contamination exists; if the flux ratio is about unity, 
indicating strong flux contamination, a large RV scattering is expected because this case is not 
considered and modeled in the current data reduction pipeline; 2, from measured spectrum, even 
the observed spectrum is a 2-D fringing spectrum in DFDI, we can still de-fringe the spectrum into 
a 1-D traditional s pectrum, on which special treatment can be performed to quantify the blending 
such as TODCOR (|Zucker et ahlboO^I ). 



7.3. Telluric Contamination 

We develop a quantitative method of estimating telluric line influence on RV uncertainty. We 
confine the discussion within the context of DFDI method, but it can be readily applied to DE 
method. We assume different levels of RV fluctuation due to atmospheric behaviors 6vn,atm (5 
m • s~^, 10 m • and 20 m • s~^). After masking a certain portion of the severely contaminated 
stellar spectrum, optimal RV uncertainty Vrms is ~3 times worse than the photon-limited RV 
uncertainty of the stellar spectrum Svrms,s- At low levels of telluric lines modeling and removing, 
it requires to mask out about 30% to 50% of a stellar spectrum to reach an optimal RV uncertainty. 
As telluric line removal becomes more effective, optimal RV uncertainty, 6vrms, reaches a smaller 
value. At a very high level of removal, in which 99% of the strength of telluric lines is removed, 
i.e., the residual of telluric line modeling is only 1%, telluric line masking becomes unnecessary. 
The entire wavelength coverage region can be used in RV measurements, and dvrms approaches the 
fundamental photon-limited RV uncertainty, Svrms,s- 
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Fig. 1. — A schematic layout of an RV instrument using the DFDI method. 
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Fig. 2. — DFDI Ilhistration. Sq{u) is a stellar spectrum; IT{v,y) is interferometer transmission; 
B{v,y) is the image taken at a 2-D detector, u is optical frequency and y is coordinate of slit 
direction. DFDI measures RV by monitoring phase shift of stellar absorption line fringes in the y 
direction (the slit direction). 
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Fig. 3. — Ex amples of synthetic M dwarf stellar spectra {V sin i=0 km • s which are generated by 



PHOENIX (jHauschildt et alJll999l : Ullard et aPboOll ). The top panel shows the spectra between 



0.8-1.35 /im, the bottom panel shows an enlarged spectral region around 1177 nm showing stellar 
line profiles. 
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Fig. 4. — Top: power spectrum of derivatives of stellar spectrum, F[dSQ/dv\, with SRF{p) = 
T[LSF{i',R)\ for different R overplotted, where LSF is line spread function; Bottom: FldSQ/du] 
shifted by = 20mm using a fixed-delay Michelson interferometer. The SRFs for different R are 
overplotted. 
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Fig. 5. — Optimal OPD correlation with Vsini (left) and spectral resolution R (right). Optimal 
OPD for /2=20,000 (solid), 50,000 (dotted), 80,000 (dashed) are used on the left panel. Vsmi=2 
(solid), 5 (dotted), 10 (dashed) km-s~^ are assumed on the right panel. Complete results of optimal 
OPD can be found in Table [TJ T^s influence on optimal OPD is not significant, rcfj=2800 K is 
adopted in the plot. 
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Resolutlon[x10^] 

Fig. 6. — Ratio of Qoptimai and Qdcviated as a function of spectral resolution, where Qoptimai is Q 
factor at optimal OPD and (^deviated is Q factor when OPD is deviated from Qoptimai by 5 mm. 
Different line styles represent different V svni while colors indicate different Tcfj. 
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Fig. 7. — Q factor as a function of spectral resolution, (left: Teg = 2400-fC; right: Teg = 3100i^. 
Open circles represent Qdfdi; crosses represent Qde! solid lines are best power-law fits for Qdfdi; 
dashed lines are best power-law fits for Qde) 
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Fig. 8. — Improvement of Qdfdi over Qde as a function of spectral resolution. 
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Fig. 9. — Comparison of QJret Qde different R. Note that Q' = Q ■ \/N~- Different color 
represents different rotational velocity. 
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Fig. 10. — Comparison of Q^fdi ^dfdi r=ioo ooo different R. Note that Q" = Q • \/N^- ■ 
The maximum of each curve is indicated by filled circle. Different color represents different 
rotational velocity, the same as Fig. [9l 
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Fig. 11. — Comparison of QJ^^ and Q'^e r=ioo ooo different R. Note that Q" = Q-\/N^- '^obj- "^^^ 
maximum of each curve is indicated by fihed circle. Different color represents different rotational 
velocity, the same as Fig. [H 
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Fig. 12. — Qdfdi as a function of V s'mi. 
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Fig. 13. — OPD optimization for IRET. OPD is optimized at R = 22, 000, V sini between 3 km • s 
and 5 km • s^^. 
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Fig. 14. — Instrument throughput 77 as a function of wavelength (800 to 1350 nm). The telescope 
transmission is 80%. It does not include atmospheric tranmission, fiber coupHng, transmission and 
seeing losses. 
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Fig. 15. — Predicted photon-limited RV uncertainty contours (in m • s~^) for IRET. The assumption 
in calculation includes: 1) texp = 30 min; 2) mj = 9; 3) ry as shown in Fig. [HI 
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Fig. 16. — S/N per pixel as a function of wavelength. The assumption in calculation includes: 1) 
^exp = 30 min; 2) mj = 9; 3) ?7 as shown in fig [Ml IRET is designed to have 4.2 pixels to sample 
one RE and 25 pixels to sample fringes in the slit direction in each spectral channel. 
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Fig. 17. — Atmosphere transmission (AT) as a function of wavelength (Generated by ATRAN ([Lord 
19921 )1. 
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Fig. 18. — Predicted RV uncertainty 5vrms as a function of the fraction of stellar spectrum used in 
RV measurements when applying the telluric line masking technique. We assume Tefj=2400 K and 
V sini of 5 km • s"^. Different levels of atmosphere RV fluctuation are indicated by different line 
styles. Svrms,s (dash-dotted line), the fundamental photon-limited RV uncertainty determined by 
an intrinsic stellar spectrum for a complete wavelength coverage from 800 to 1350 nm, is plotted 
for comparison. 
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Fig. 19. — Predicted RV uncertainty as a function of fraction of stellar spectrum used in RV mea- 
surements when applying the telluric line modeling and removing technique. We assume Tgff =2400 
K and V sini of 5 km • s~^. Asterisks connected by solid line represent the fundamental photon- 
limited RV uncertainty. Different line styles represent different level of telluric lines removal (i.e, 
removed telluric line strength). 6vn,atm=W m • is assumed in the calculations. 
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Fig. 20. — Predicted optimal RV uncertainty as a function of removed telluric line strength. Dashed 
line represents the fundamental photon-limited RV uncertainty, Vrms,s- 
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Table 1. OPD choice as a function of R and V sini at different T^s 



R 





1 




2 




3 


4 


V 


sin i 
5 


[km 


■ s 
6 


-'1 


7 


8 


9 


10 


5,000 


19,27,29 


19,19,28 


19 


19 


19 


15,15,17 


15,15,15 


13 


13,13 


11 


12 


12 


10,11,11 


10,10,10 


10.10,10 


10,10,10 


10,000 


21,23,27 


19,21,26 


19 


19 


20 


17,17,17 


15,15,15 


13 


13,13 


11 


12 


12 


11,11,11 


10,10,10 


10.10,10 


10,10,10 


15,000 


21,23,26 


21,22,23 


19 


20 


20 


17,17,17 


15,15,15 


14 


14,14 


13 


13 


13 


12,11,11 


11,11,11 


11.10.10 


10,10,10 


20,000 


22,23,26 


21,23,24 


20 


21 


21 


19,19,19 


17,17,17 


15 


15,15 


14 


14 


14 


13,13,13 


12,12,12 


11,11,11 


11,11,11 


25,000 


23,24,26 


23,24,25 


21 


22 


22 


19,20,20 


19,18,18 


17 


17,16 


16 


15 


15 


14,14,14 


14,13,13 


14.13.12 


12,12,12 


30,000 


24,26,26 


24,25,26 


23 


23 


24 


21,21,21 


19,19,19 


18 


17,17 


17 


16 


16 


16,15,15 


14,14,14 


14,14.14 


14,14,14 


35,000 


26,26,28 


26,26,26 


24 


24 


25 


22,22,23 


21,21,20 


19 


19,19 


19 


18 


17 


17,17,17 


16,16,16 


16,16.16 


16,16,15 


40,000 


27,28,28 


27,27,28 


26 


26 


26 


24,24,24 


22,22,22 


21 


20,20 


19 


19 


19 


19,18,18 


19,18,18 


19,18,18 


19,18,16 


45,000 


29,29,30 


28,28,29 


27 


27 


28 


25,25,26 


24,23,23 


22 


22,21 


22 


20 


20 


22,20,20 


22,20,18 


22,18,18 


22,18,18 


50,000 


30,30,31 


29,30,30 


28 


28 


28 


27,26,26 


24,24,24 


24 


22,23 


22 


22 


22 


22,22,22 


22,22,22 


22,22.22 


24,22,22 


55,000 


32,32,32 


31,31,32 


30 


30 


30 


27,28,28 


27,26,26 


24 


24,24 


24 


24 


24 


24,24,22 


24,24,22 


24.24.22 


24,24,24 


60,000 


32,32,34 


32,32,33 


32 


31 


32 


29,28,28 


27,27,26 


27 


26,26 


24 


24 


24 


24,24,24 


24,24,24 


24.24.24 


24,24,24 


65,000 


34,34,35 


34,34,34 


32 


32 


32 


32,30,30 


29,28,28 


27 


27,26 


27 


26 


26 


27,24,26 


24,24,26 


24,24,24 


24,24,24 


70,000 


35,35,36 


35,35,36 


34 


32 


34 


32,32,32 


32,30,30 


32 


28,28 


27 


28 


28 


27,28,26 


27,24,26 


24,24,26 


24,24,26 


75,000 


37,36,37 


37,36,37 


35 


35 


35 


32,32,32 


32,32,32 


32 


32,32 


32 


32 


28 


32,32,28 


32,32,32 


32,32,32 


32,24,32 


80,000 


37,37,39 


37,37,37 


37 


36 


36 


35,34,35 


35,32,32 


32 


32,32 


32 


32 


32 


32,32,32 


35,32,32 


35,32,32 


35,32,32 


85,000 


40,39,41 


40,39,39 


37 


37 


37 


37,36,36 


35,32,32 


35 


32,32 


35 


32 


32 


35,32,32 


35,32,32 


35,32,32 


35,32,32 


90,000 


40,41,41 


40,41,41 


40 


37 


39 


37,36,36 


37,36,36 


37 


36,36 


35 


32 


36 


35,32,36 


37,37,36 


37,37,36 


37,37,36 


95,000 


41,41,41 


41,41,41 


40 


41 


41 


40,37,37 


37,37,36 


37 


37,36 


37 


37 


36 


37,37,36 


37,37,37 


37,37,37 


37,37,37 


100,000 


41,41,41 


41,41,41 


40 


41 


41 


40,40,41 


37,37,37 


37 


37,36 


37 


37 


36 


37,37,37 


37,37,37 


37,37,37 


37,37,37 


105,000 


41,41,41 


41,41,41 


41 


41 


41 


40,41,41 


40,40,41 


37 


37,37 


37 


37 


37 


37,37,37 


37,37,37 


37,37,37 


37,37,37 


110,000 


41,41,41 


41,41,41 


40 


41 


41 


40,41,41 


40,41,41 


37 


41,41 


37 


37 


41 


37,37,37 


37,37,37 


37,37,37 


37,37,37 


115,000 


41,41,41 


41,41,41 


41 


41 


41 


40,41,41 


40,41,41 


37 


41,41 


37 


41 


41 


37,37,41 


37,37,37 


37,37,37 


37,37,37 


120,000 


41,41,41 


41,41,41 


41 


41 


41 


40,41,41 


40,41,41 


37 


41,41 


37 


41 


41 


37,37,41 


37,37,37 


37,37,37 


37,37,37 


125,000 


41,41,41 


41,41,41 


41 


41 


41 


40,41,41 


40,41,41 


37 


41,41 


37 


41 


41 


37,37,41 


37,37,37 


37,37,37 


37,37,37 


130,000 


41,41,41 


41,41,41 


40 


41 


41 


40,41,41 


40,41,41 


37 


41,41 


37 


41 


41 


37,37,41 


37,37,37 


37,37,37 


37,37,37 


135,000 


41,41,41 


41,41,41 


40 


41 


41 


40,41,41 


40,41,41 


37 


41,41 


37 


41 


41 


37,37,41 


37,37,37 


37,37,37 


37,37,37 


140,000 


41,41,41 


41,41,41 


40 


41 


41 


40,41,41 


40,41,41 


37 


41,41 


37 


41 


41 


37,41,41 


37,37,41 


37,37,37 


37,37,37 


145,000 


41,41,41 


41,41,41 


40 


41 


41 


40,41,41 


40,41,41 


40 


41,41 


37 


41 


41 


37,41,41 


37,41,41 


37,41,41 


37,41,41 


150,000 


41,41,41 


41,41,41 


40 


41 


41 


40,41,41 


40,41,41 


40 


41,41 


37 


41 


41 


37,41,41 


37,41,41 


37,41,41 


37,41,41 



Note. — The first number in each tab is the optimal OPD for T^ff — 2400/C, the second number is for T^^f — 2S00-fC and the third one is for 
Tjjff — 3100-fC. OPD is in the unit of mm. OPD ranging from 10mm to 41mm is considered in calculation. 



Table 2: Power Law Index x as a function of Spectral Resolution R for DFDI and DE (T^ff = 2400K) 









DFDI 






DE 




V sini 
[km • s^"'"] 


R 


5,000-20,000 


20,000-50,000 


50,000-150,000 


5,000-20,000 


20,000-50,000 


50,000-150,000 







0.62 


0.59 


0.31 


1.08 


0.93 


0.44 


2 


X 


0.63 


0.56 


0.27 


1.07 


0.89 


0.38 


5 


0.62 


0.45 


0.16 


1.01 


0.69 


0.21 


10 




0.58 


0.28 


0.09 


0.87 


0.39 


0.10 
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Table 3. Spectral Resolution and wavelength coverage on a given 

detector 



R 


AA 


Amin Amaj^ 




(nm) 


(nm) 


25,000 


480 


800-1280 


30,000 


400 


800-1200 


40,000 


300 


850-1150 


50,000 


240 


880-1120 


60,000 


200 


900-1110 


70,000 


170 


910-1080 


80,000 


150 


920-1070 
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Table 4: Q" comparison of DFDI and DE as a function of V sin i 

DFDI DE 





V sa\i [km • s "'"] 


Roptinial 


Qdfdi 


Roptimal 




Qdfdi/ Ode 







50,000 


7502 


110,000 


6623 


1.065 


a=0.5 


2 


50,000 


6806 


75,000 


6001 


1.134 










5 


25,000 


4979 


50,000 


4424 


1.125 




10 


15,000 


3394 


25,000 


2996 


1.133 


a=1.0 



2 
5 

10 


5,000 
5,000 
5,000 
5,000 


26384 
24297 
18929 
12877 


30,000 
25,000 
15,000 
10,000 


8705 
8490 
7884 
7239 


3.031 
2.862 
2.401 
1.779 



Table 5. Predicted IRET Performance 



Name 


nij 


Toff 




K 


Sv a'' 


Knz" 






(K) 


(km-s"i) 


(m-s-i) 


(m-s-1) 


(m-s-i) 


GJ 1214 


9.75 


3000 


2 


12 


2.4 


1.0 


GJ 176 


6.46 


3500 


1 


4.1 


0.58 


0.57 


GJ 179 b=* 


7.81 


3400 


1 


26 


1.05 


0.66 


GJ 436 b-* 


6.9 


3684'' 


1 


18.7 


0.71 


0.59 


HIP 57050 b^ 


7.61 


3190 


1 


38 


0.91 


0.68 


GJ 649 b'' 


6.45 


3700^* 


1 


12 


0.58 


0.54 



Note. — a: V sini is assumed to be 1 km ■ s^^ if otherwise specified in references; 
b: the fundamental photon-limited RV uncertainty; c: velocity semi-amplitude if 
there is a habitable Earth-like planet locating at 0.05 AU from a host star; d: we 
assume Tofi to be 3500K because we do not have synthetic stellar spectrum with T^a 
higher than 3500K. 
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Table 6. Mid-Late M Dwarfs Available for IRET" 



mj 


Number of Targets'" 


J<6 


25 


J<7 


97 


J<8 


270 


J<9 


868 


J<10 


2458 


J<11 


6267 


J<12 


12775 







Note. — a; LSPM Catalog (|Lepme fc Shara 
cut (V-J>3) and mj cut. 



20051 ) : b: we apply both color 



